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ABSTRACT 

It has been recently shown that the relative velocity between the dark matter and the baryons 
^bc) at the time of recom bination can affect the structure formation in the early universe 
Tseliakhovich fc Hiratal l2010f) . We statistically quantify this effect using large cosmological simu- 
lations. We use three different high resolution sets of simulations (with separate transfer functions for 
baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity 
between dark matter and baryons. We show that the total number density of halos is suppressed 
by ~ 20% at z = 25 for v^ c — lovbc, where cr v bc is the variance of the relative velocity, while for 
^bc = 3.4cr v bc the relative suppression at the same redshift reaches 50%, remaining at or above the 
30% level all the way to z = 11. We also find high abundance of_"empty halos", i.e., halos that 
have gas fraction below half of the cosmic mean baryonic fraction /b. Specifically we find that for 
fbc = lcVbc oil halos below 10 5 M Q are empty at z > 19. The high abundance of empty halos results 
in significant delay in the formation of gas rich mini-halos and the first galaxies. 



1. INTRODUCTION 

Observations of the cosmic microwave background 
(CMB) show that the Universe at cosmic recombination 
(redshift z ~ 10 3 ) was remarkably uniform apart from 
small spatial fluctuations in t he energy density and in the 
gravitational potential (e.g.. [Bennett et"aTI[l9 96). The 
primordial inhomogeneities in the density distribution 
grew over time and eventually led to formation of galax- 
ies, galaxy clusters and large-scale structure. Various 
physical processes contributed to the gro wth of pertur- 
bations in dark matter and ba ryons (e.g., lPeebleslll9861 . 
119931 : iMa fc Bertschingedll995l ). which evolved very dif- 
ferently before recombination. As the result, baryon den- 
sity fluctuations were much smaller than fluctuations in 
the dark m atter density at the time of rec ombination. 

Recently Tseliakhovic h fc Hiratal ((2010) showed that 
not only the amplitudes of the dark matter and baryon 
density fluctuations were different but also were their ve- 
locities. After recombination, the sound speed of baryons 
dropped down to thermal velocities, while the dark mat- 
ter velocity remained unaffected - thus, the relative ve- 
locity of the dark matter with respect to the ba ryons 
become supersonic. iTseliakhovich fc Hiratal (|2010f ) also 
show that the relative velocity between the dark matter 
and baryons is coherent on scales of a few mega-parsec 
and is in the order of ~ 30 km sec^ 1 at the time of cosmic 
recombination. This relative velocity is often called the 
"stream velocity" in the literature, and throughout this 
paper we will use this term. The stream velocity effect 
was previously overlooked, because the velocity terms are 
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actually of the second order in perturbation theory. The 
stream velocity effect leads to power suppression at mass 
scales that correspond to masses of the first bo und ob- 
jects in the Univers e fe.g.. | Yoshida et al.ll2003a[ ). Using 
the Press-Schechter dPress fc Schechtedll974T ) formalism, 
ITseliakhovich fc Hiratal (|2010D showed that the number 
density of halos is suppressed by more than 60% for ha- 
los with M = 10 6 Mo a t z = 40. In a subsequent pa- 
per, ITseliakhovich et ail (120101) also included th e bary- 
onic component following [Naozl&^arjfijiia] ((lOOU). They 
found that the stream velocity also results in much higher 
minimum mass for a gas rich halo at high redshifts, as 
compared to the case without the stream velocity (e.g. 
iNaoz fc Barkmall2l)u7l). 

ITseliakhovich fc Hiratal ()2010D and 

ITseliakhovich et al.l (J2010) findings have prompted 
several studi es on the stream ve locity effect in numerical 
simulations. iStacv et al.l (|2011l ) studied the collapse of 
gas in minihalos and found only a mild delay (~ 10 7 yr) 
for the time of collapse at low redshift. They also 
estimated the difference in redshifts of collapse for the 
minimum halos that can collapse, with and without the 
stream velocity effects, and found only minimal effect 
due to the thermal evolution of the gas collapsing into a 
minihalo at low redshifts. Their results also show that 
the delay of the struct ure formation beco mes significant 
with higher redshift. Mai o et "al~l (j 2 1 lh examined the 
influence of the stream velocity on the formation of 
the smallest bound objects. They showed that the 
abundance of mini-halos (< 10 7 M Q ) is reduced by a 
few percent and that the baryon fraction in those halos 
is red uced by somewhat larger amount. iDalal et all 
(2010) showed that the characteristic imprint of the 
stream velocity effect on the power spectrum of the 
earliest structures can be used to distinguish the effects 
of minihalos on intergal a ctic g as before and during 
reionization. iGreif et all (|2011| ) used high resolution, 
moving-mesh simulations to study the effect of different 
initial stream velocities on the virialisation of the gas 



2 



Naoz et al. 



in minihalos; their simulations suggest that the stream 
velocity substantially delay the onset of gravitational 
collapse in minihalos. Although all these studies 
consistently show some delay in the collapse of early 
gas clouds, a thorough statistical calculation is clearly 
needed to quantify these results, and, especially, to 
examine critically the effect on early baryonic structure 
formation. 

In this paper we quantify the suppression of structure 
formation due to the stream velocity effect by means of 
cosmological simulations. We use three different sets of 
high resolution simulations in order to study the stream 
velocity effect systematically, thus understanding the 
overall picture (instead of co ncentrating on sp ecific ha- 
los). In a companion paper (|Naoz et al.ll201ll in prep.) 
we explore the effect of stream velocity on the gas fraction 
in halos and compare simul ation results to the predic - 
tions from linear theory fe.g. iTseliakhovich et al.ll2010h . 
We use a set of simulations with different box sizes, par- 
ticle numbers, and the values for the stream velocity to 
analyze the suppression of the structure formation as a 
function of the stream velocity. 

We first describe the parameters and initial conditions 
of our simulations ($2]). We present our results and anal- 
ysis of the halo mass functions in the simulations in Sec- 
tion [3] Finally we offer a discussion in gl 

Throughout this paper, we adopt the following cosmo- 
logical parameters: (J7a, ^m, O;,, n, <7 8 , Hq)= (0.72, 0.28 , 
0.046, 1, 0.82, 70 km s~ 4 Mpc^ 1 ) (jKomatsu et al.ll2009h . 

2. COSMOLOGICAL SIMULATIONS 

2.1. Simulation Setup 

We use parallel AT-body /hydrodynamics solver 
GADGET-2 (jSpringel et al.1 12001 ISprineell [2005h . We 
describe the general features of the simulations, which 
are also summarized in table [TJ 

1. The first set, named U N = 256", uses a total of 
2 x 256 3 dark matter and gas particles within a 
cubic box of 200 comoving kpc on a side. To real- 
ize statistically significant number of halos in such 
a small box, we artificially increase gravitational 
clustering in the simulation by setting eg — 1.4. 
We choose this box size so that a 10 4 M halo is 
resolved with ~ 500 particles - the valu e needed to 
estim ate the halo gas fraction reliably (jNaoz et al.l 
2009). The gravitational softening is set to be 
40 comoving pc, well below the virial radius of a 
10 4 Mq halo (~ 680 comoving pc). All the simula- 
tions in this set are initialized at z = 199. 

2. The second set, named "A = 512", uses a total 
of 2 x 512 3 dark matter and gas particles within 
a cubic box with the size of 700 kpc. In this set 
we also artificially increase <r$ to 1.4. With these 
parameters, a halo with 500 dark matter particles 
has a mass of - 5 x 10 4 M Q . The softening length 
is set to be 68 comoving pc. All the simulations in 
this set are initialized at z = 199. 

3. The final set of simulations uses 2 x 768 3 dark mat- 
ter and gas particles (which we name the "A = 
768" set) in a 2 Mpc box, and starts at z = 99. 
For these parameters a halo with 500 dark matter 



TABLE 1 
Parameters of the simulations 



SIM 


l>bc,0 


<r v bc stream 


Shift in 




km sec 


velocity 


position 




256 runs, 


0.2 Mpc, 


^in — 199 


256o 








No 


256io-+p 


5.8 


1 


Yes +18.05 kpc/h 


256i. 7ct 


10 


1.7 


No 


256i.7 CT+p 


10 


1.7 


Yes +30.7 kpc/h 


2562.6<t 


15 


2.6 


No 


2563.4<T 


20 


3.4 


No 




512 runs, 


0.7 Mpc, 


%vn — 199 


512 








No 


512 iCT 


5.8 


1 


No 


512x.ro. 


10 


1.7 


No 


5123.4,7 


20 


3.4 


No 




768 runs, 


2 Mpc, 


^iti — 99 


768 








No 


768^ 


3 


1 


No 


768i.7<7 


5 


1.7 


No 


768 3 .4<7 


10 


3.4 


No 



particle has a mass of - 10 5 M Q . The softening 
length is set to be 0.2 comoving kpc. We use the 
"correct" value of <j$ — 0.82 for this simulation set. 



In each simulation set, we explore a variety of the values 
for the stream velocity (see £12.31 and table [T]). 

2.2. Initial Conditions 

As has been show n by iNaoz et"ai1 ([2006?) and 
iNaoz fc Barkanal (|2007ft . setting up initial conditions on 
the small spatial scales is a delicate issue. High accu- 
racy in initial conditions is crucial for accurately pre- 
dicting the halo mass function in the lowest mass regime 
(M < 10 7 M p ). 

Following INaoz et"aLl ()2010[ ). we generate separate 
transfer fu nctions for dark matter and baryons as de- 
scribed in INaoz fc Barkana (|2005l ). Ideally, streaming 
velocities should be realized in the initial conditions in 
a self-consistent way with the transfer functions that are 
calculated with streaming effect. However, these trans- 
fer functions do not take into account the stream veloc- 
ity effect, in consistency with the previous calculations 
that included the stream velocity in simulations. For all 
runs, glass-like initial conditions were generated using 
the Zel'dovich approximation. For baryons, we have used 
a glass file with positions shifted by a random vector, 
thus removing artificial coupling between ne arby dark 
matter and gas particles ()Yoshida et a l. 2003b). We note 
that we have used the same phases for dark matter and 
baryons in all of our simulations except for those with 
position shift, see §2.4| . 

2.3. Stream Velocity in the Simulations 

iTseli akhovicb Tfc" Hiratal ([2010D showed that, while the 
stream velocity varies in space, its coherence length is 
quite large, many mega-parsecs. Hence, on scale of our 
simulation boxes it can be treated as constant bulk mo- 
tion of baryons with respect to dark matter. We include 
the effect of stream velocity by adding, at the initial red- 
shift, an additional velocity to the x component of the 
baryons velocity vector. We test a range of values for the 
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stream velocity, which is convenient to quantify in terms 
of its rms value on small scales, er v bc- Specifically, we test 
«bc = lcvbc through Vhc = 3.4cr v b c for all the simulations 
sets (see table Q]). 

2.4. Position Shift in the Initial Conditions 

It is customary to set the initial conditions such that 
the phases of dark matter and baryons are the same, 
as dictated by the linear perturbation theory. However, 
because the stream velocity effect is of the second order, 
it violates the assumption. Namely, since the baryons 
move in the dark matter reference frame, the respective 
density perturbations that were co-located at the time 
of recombination, become separated in space. In other 
words, the phases of dark matter and baryons at the 
same spatial location on sufficiently small scales become 
unrelated. 

The relative shift between dark matter and baryons is 
easily calculated through the geodesic equation, 



100000 



ax 



Vbc(a) 



(1) 



where x is the comoving distance, a is the scale factor, 
and ^bc(a) is the lagrangian velocity of a baryonic fluid 
element in the dark matter reference frame as a function 
of the scale factor a. Hence 



. , f* Vb c (a) 
Axbc{a) — I at 



(2) 



where t is the cosmic time at the scale factor a. Tak- 
ing this integral, we find the relative shift in the posi- 
tion between baryons and dark matter at z- m = 199 is 
Axhc = 18.5 comoving kpc/h for vt> c = loVbc and is 
~ 30.7 comoving kpc/h for v hc = 1.7a v b c - We test the 
effect of the position shift on two of the TV = 256 runs 
by adding the above values to the x component of the 
baryons at the initial redshift z- ln , thus effectively chang- 
ing the phases of the baryons relative to dark matter. 

As we find below, the relative shift between baryons 
and dark matter is a small effect, and it does not invali- 
date the rest of our simulations, where the position shift 
is neglected. 

2.5. Halo definition 

We locate dark matter halos by running a friends-of- 
friends group-finder algorithm with a linking parameter 
of 0.2 (only for the dark matter component). We use 
the identified particle groups to find the center of mass 
of each halo. After the center is located, we calculate 
density profiles of dark matter and baryons separately, 
assuming a spherical halo and using 2000 radial bins be- 
tween r m i n = kpc and r max = 20 kpc. Using the density 
profiles, we find the virial radius r V i r at which the total 
overdensity is 200 times the mean background density, 
and compute the mass and the gas fraction of each halo 
within that radius. 

Recently, IMore et al.l (|2011l ) showed that halos identi- 
fied by the friends-of-friends algorithm enclose an average 
overdensity that is substantially larger than 200 and its 
specific value depends on the halo concentration. In our 
approach we use the friends-of-friends algorithm only for 
finding the center of mass of a halo, and compute the 
actual halo mass using the spherical overdensity of 200. 




Fig. 1. — Cumulative halo mass function for different simulations 
and valus of Vb c at z = 15. We show the N = 256 (long dashed 
lines), N = 512 (solid lines), and N = 768 (short dashed lines). We 
consider the values for the stream velocity of V), c = 0, i>b c = lCvbo 
v hc = 1.7cr v b c , and v hc = 3.4<r vbc , from top to bottom in each sets 
of lines. 

We only retain halos that contain at least 500 dark 
matter particles within their virial radii. The choice al- 
lows us to estimate halo masses to about 15% precision 
(jTrenti et al.l 12010(1 and to estimate halo gas fraction s 
reliably to a similar level of accuracy ([Naoz et al.ll2009h . 

3. RESULTS 
3.1. Suppression of the Halo Mass Function 

We first consider all the simulations that do not in- 
clude the relative position shift between baryons and 
dark matter (i.e., all runs except N = 256_l £ t+ p and 
N = 256i. 7<7+p ). iTseliakhovich fc Hiratal l|2010j ) calcu- 
lated the number densities of collapsed halos with and 
without the stream velocity, and illustrated the stream 
velocity effect on the abundance of small halos. Specif- 
ically, based on the Press-Schechter formalism, they 
showed that the number density of haloes with the rms 
stream velocity A v b c (> M) is suppressed by more than 
60% at the mass scale of M ~ 10 6 Mg in comparison to 
the case with no stream velocity (N (> M)). 

In Figure Q] we show the halo mass functions at z = 
15 for all our simulations. Indeed, as can be seen, the 
halo number density is a strong function of the stream 
velocity, resulting in substantial suppression of the halo 
mass functi on for typical values of t h e stre am velocity^ 

Following Tseliakhovi ch fc Hiratal (|2010[ ). we quantify 
the suppression of the halo mass function as 



_ N vhc {> M) - N {> M) 
N (> M) 



(3) 



A more relevant quantity is, perhaps, the abundance of 
halos that contain substantial amount of gas. It is such 
halos that may host first stars and first supernovae, and 
also serve as sinks of ionizing radiation during reioniza- 

6 The lower halo mass function for the N = 768 runs is the result 
of our simulation setup; N = 768 have much lower <rg than other 
simulation sets. 
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Fig. 2. — Effect of the stream velocity on the cumulative halo 
mass function for the N = 768 runs, for redshifts of (from bottom 
to top panels) 19, 15, and 11. The relative difference for all halos 
(A v ) is shown with the solid lines, while dashed lines show the 
relative difference for non-empty halos (A v j g ), for various values of 
the stream velocity v^ c = l<r v b c , fbc = 1 -f&vbc, an d ^bc = 3.4cr v b c 
(purple, brown and red lines respectively). Note that A v j for 
■"be = 3.4cr v b c at z = 19 is exactly —1 (i.e., there are no halos with 
gas fraction above /b/2 in that simulation). The opposite is true 
for Dbc = lo" v bc at z = 15 and z = 11 - there are no empty halos in 
those runs. The dotted vertical line in each panel marks the mass 
at which the halo number density is 5 kpc~ . 

tion. Therefore, we also introduce 

A _ iVvbc(> M, f g > / b /2) - JV (> M, f g > / b /2) 
Vjg N (> M,f g > f h /2) 

(4) 

as a relative difference between the mass functions of 
halos that contain more gas_ than one half of the cos- 
mic mean baryon fraction (/b) in simulations with and 
without the stream velocity, iV v b c (> M,f g > /b/2) and 
N (> M, f g > /b/2) respectively. 

We would call halos with gas fractions less then /b/2 
"empty halos", although, admittedly, this definition is 
somewhat arbitrary and a different gas fraction threshold 
could also be considered. 

First we show in Figure [5] the relative suppression 
of halo mass function due to the stream velocity for 
all and for non-empty halos only (A v and A Vi / ; solid 
and dashed curves respectively) for the N = 768 set 
of simulations. At z = 19 the suppression of the total 
number density of halos is about 50% for 10 6 M for 
u bc = 3.4er vbc and ~ 15% for Vh c — lfvbc- As time 
goes by, the effect of the stream velocity diminishes, but 
remains clearly visible even at the lowest redshift we con- 
sider, z = 11, still present. 

Perhaps the most interesting feature shown in Figure 
[5] is the high abundance of empty halos. At z = 19, 
more than 20% of the halos below 10 6 M are empty 
for «bc = l-7<r v bc and all of the halos in the simula- 
tion are empty for u bc = 3.4er v b c . These empty halos 
lay above the minimum cooling mastO] (~ 3 x 10 5 Mq, 

7 We note that we in the simulations presented in the paper we 
included only adiabatic cooling. 



Fig. 3. — Same as Figure |21 for the TV = 512 runs at redshifts 
(from bottom to top panels) 25, 19, and 15. 
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Fig. 4. — Same as Figure [2] for the TV = 256 runs at redshifts 
(from bottom to top panels) 25, 19, and 15 and for 4 values of the 
stream velocity v^ c = lovbci fbc = l- 7(T vbci ^bc = 2.6<r v b c , and 
v bc — 3.4cr v b c (purple, brown, green and red lines respectively). 

e.g.. lTrenti fc Stiavellil[2009]) . Thus, we conclude that, in 
patches of the universe where the stream velocity is high, 
the formation of the first stars and galaxies is delayed. 
This behavior is consistent with the N = 512 simulation 
set (Figure [3j) , where, for t> bc = 3.4er v b c , most of the ha- 
los below ~ 10 6 Mq are empty of baryons at z = 19, and 
all of the halos in the simulation are empty at z = 25. 
Furthermore, for «bc = l-7cr vbc , more than 50% of the 
halos below ~4x 10 5 M Q are empty for z > 19. 

3.2. Saturation of the Halo Mass Function Suppression 
in the Low Mass Limit 

As expected, larger values for the stream velocity re- 
sult in stronger suppression of the halo mass function. 
However, in the N = 256 simulation set the effect clearly 
saturates for high values of Ub c (see Figure |4|. In order 
to explore the nature of this saturation further, we com- 
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Fig. 5. — Halo mass function suppression factor A v at z = 15 for 
all our simulation sets for different values of the stream velocity 
(from bottom to top) of v^, c = l<r v bcj 1.7<7 v b c , and 3.4(r v b c . In 
all panels N = 768, N = 512, and N = 256 simulation sets are 
represented by dashed, solid, and dotted lines respectively. All 3 
simulations sets are reasonably consistent with each other within 
the cosmic variance. Note that to check consistency with somewhat 
lower mass values, we have lowered the threshold for the number 
of particles per halo to 300 particles per halos only for the N = 768 
runs and only in order to generate this figure. 

pare A v at z = 15 (the latest epoch that our N = 256 
runs can be continued to) between different simulation 
sets in Figure [SJ 

Overall, the agreement between all three simulation 
sets is reasonable! N = 256 set suffers from poor statis- 
tics for M > 10 5 M Q , so the agreement between that se- 
ries and two other sets is somewhat worse. Nevertheless, 
the effect of saturation for low masses is clearly visible 
for all values of Ubc- 

The reason why A v as a function of halo mass is ex- 
pected to saturate at the low mass limit is understood 
as follows. For halo mass of 10 5 M Q the escape veloc- 
ity is about 0.77 km sec -1 , while the stream velocity for 
Ubc = 3.4cr v bc at z = 15 is 1.6 km sec -1 and at z — 25 is 
2.6 km sec -1 . Thus, it is not surprising that halos below 
10 5 Mq are empty in that redshift range - the stream 
velocity is simply much larger then the halo escape ve- 
locity, so the dark matter halo is unable to accrete any 
gas. Hence, such a halo will have only about (1 — /b) of 
the mass it would have had if v^c = 0. Since the cumula- 
tive mass function N(> M) oc Af -1 in this mass regime, 
the loss of baryons by low mass halos results in /b « 0.17 
reduction in the halo mass function. 

At larger masses the interaction between the baryons 
moving with the highly supersonic velocity and the dark 
matter halo becomes more complex. For example, in 
the Ubc = 3.4cr v bc case the mass function is smaller by 
~ 40% at 1O 5 M < M < 3 x 10 6 M Q . This magnitude 
of the effect cannot be explained by the loss of baryons 
alone; some non-negligible fraction of dark matter has to 
be lost by a halo as well. Note that the simulations sets 
N = 512 and N — 768 have different erg, but agree well 

8 Note that this agreement is present through all the simulations 
despite the different ag assumed. 



Fig. 6. — Halo mass function suppression factor A v as a function 
of time for 5 X 10 s Mq halos. In this figure we only show N = 768 
(dot-dashed curves) and N = 512 (solid curves) simulation sets. 
We also show the suppression factor for non-empty halos, A v f , 
with dotted lines. We consider i>b c = lo" v b c , "be = l-7 CT vbc> an d 
"be = 3.4<r vbc from bottom to top panels respectively. 

in the magnitude of the halo mass function suppression. 
Hence, these results suggest that the suppression effect 
is largely due to non-trivial non-linear interaction of the 
baryonic flow and the dark matter halo. 

In Figure [5] we consider the time evolution of the halo 
mass function suppression factor for a halo mass that 
lays safely above the molecu lar hydrogen cooling mass 
(e.g.. iTrenti fc Stiavelli 2009|). We choose a value of 5 x 
10 5 M , which is well resolved and statistically reliably 
represented in both N = 768 and N = 512 simulation 
sets. This comparison suggests significant delay in the 
formation epoch of gas rich halos, which, in principle, 
could have formed stars. We find that for v^ c = lavbc 
the sup pression is modes t , in agreement w ith previous 
studies (|Maio et al.|[20Tlt iStacv et ail 1201 lft . However, 
for larger values of the stream velocity (vbc = 1.7<r v bc 
and Vbc = 3.4er v bc) we find a significant delay and large 
abundance of empty halos. 

3.3. Position Shift 



As we discussed in §2.41 the stream velocity also causes 
the baryons to shift in position relative to dark mater. 
We have tested this effect in two of the N = 256 simula- 
tion runs. In Figure [7] we show the relative suppression 
of halo formation comparing between the cases with and 
without the position shift (thin and thick curves respec- 
tively). By z — 15 the relative shift in the positions is 
erased; at earlier times the effect is noticeable, but it 
never becomes the dominant effect in the change in the 
halo mass function induced by the stream velocity. 

3.4. Suppression of the Clumping Factor 

The baryonic clumping factor plays an important 
role in the penetration and escape of radiation from 
an inhomogeneous medium and thus the effects of the 
stream velocity also depicted in the suppre ssion of the 
clumping factor (e.g.. iBarkana fc Loe b 2001). Adopting 
iSpringel &: Hernquistl ([20031 ) for estimating the clumping 
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Fig. 7. — Halo mass function suppression factors A v (solid 
curves) and A Vj * (dashed curves) for t) bc = l<x v bc and v bc = 
1.7er vbc cases (purple and brown curves respectively). Thin curves 
represent simulations that properly account for the position shift 
(256i CT +p and 256i.7 CT +p runs), while thick curves show our fidu- 
cial approach, without adding a position shift to the initial con- 
ditions (i.e., initializing the simulation with the same phases for 
the baryons and dark matter, 256i CTvbc and 256i.7 CT ). As one can 
see, the relative shift in the positions of baryons and dark matter 
makes only a modest correction. 



TABLE 2 

Fitting Parameters 



"be 


ai 


a.2 


«3 


«4 





415.1 ±66.4 


-0.32 ± 0.01 


0.06 ±0.12 


0.13 ±0.09 




393.8 ±42.3 


-0.33 ±0.01 


0.05 ± 0.07 


0.15 ± 0.08 


1.7<T vbc 


439.4 ±54.4 


-0.35 ± 0.01 


0.12 ± 0.15 


0.1 ±0.06 


3.4<r vbc 


607.1 ± 75.7 


-0.41 ± 0.01 


0.55 ± 0.27 


0.04 ± 0.02 



factor C from SPH simulations we define 

„ Ei m MM< &Pc)~ 1 Y JO m b.jpA< A Pc) 

where nib is the mass of the baryon particles in our sim- 
ulations, p c is the critical density in the Universe, and 
the summat ion is over all gas particles. 

Following iMaio et all (2006) we sum over all gas par- 
ticles whose density pi is smaller than a given thresh- 
old Ap c . Thus, we avoid numerical artifacts that are 
produced by gas particles belonging to collapsed objects 
which would artificially increase the value of C due to 
their high density. In Figure [8] we show the clumping 
factor for the TV = 768 run. We choose to show the re- 
sults of only this set of runs because the clumping factor 
depends on the assumed initial fluctuation amplitudes. 
For this run the gas particle mass is nib = 225.5 Mlq for 
all SPH particlefl In Figure bottom panel, we show 
the clumping factor for overdensity of A = 200, and we 
also test A = 500 (dotted lines). Note, that we omit the 

9 Unlike the variables A„ and A„ j the clumping factor is very 
sensitive to the assumed initial fluctuation amplitudes (because 
essentially we count the number of clumps, which is larger for large 
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Fig. 8. — The clumping factor as a function of redshift. In the 
bottom panel we consider the clumping factor evolution for v bc = 
0, v hc = lcr vbc , v bc = 1.7o- v bc> and v bc = 3.4o- vbc (blue, purple, 
brown and red lines respectively), choosing A = 200. We show 
in dashed lines the resulted fit to equation J|>]| ; see table [2] for the 
fit parameters. In the top panel we consider the suppression ratio 
C vbc /C vbc — o — 1 (see text) for overdensity case of A = 200 (we use 
the same color code as in previous figures). Note that the results 
for z = 13 and z = 10 of the run 768o are missing due to computer 
failure. 

A = 100 lines from the figure, to reduce the clutter, this 
choice gives a typical 1% reduction of the clumping fac- 
tor at z = 14 and about 15% at z = 10 (for all v b c cases). 
We also consider (top panel) the relative suppression ra- 
tio C v bc/C v bc=o — 1) where C v bc is the clumping factor 
for the case with stream velocity and C v bc=o is for the 
case without the stream velocity effect. 

The growth of the clumping factor as a function of time 
is quite regular with a behavior, at high redshift, that 
resembles exponential asymptotic decay. We assume an 
exponential behavior and find a simple fit in the form of: 

C(z) = a ie a2Z + a 3 exp a4Z . (6) 

We summarize the fit results (for A = 200 for different 
values of t>b c ) in table [2] 

As shown in Figure [8] the stream velocity suppresses 
the clumping of baryon especially at lower redshift. At 
high redshift the clumping is similar and indeed low in 
all cases. However, as time goes by, the dumpiness in 
the case without stream velocity increases dramatically, 
while it is smaller by about 50% for v b c — 3.4cr v b c - Figure 
[8] also shows that the suppression reaches a saturation 
at z — 14 — 12 and below this redshift we find smaller 
differences. 

4. DISCUSSION 

We have used three-dimensional hydrodynamical simu- 
lations to investigate the effects of stream velocity on halo 
mas s function in the early universe (in a companion pa- 
per iNaoz et al"1l2011l we study the effect on the g a s frac - 
tion and filtering mass). [Tseliakhovi ch fc Hiratal (|20100 
showed recently within the frame work of linear theory, 
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that the initial velocity difference between baryons and 
dark matter after recombination significantly suppresses 
the halo mass function at small scales. 

We found that the total halo mass function as a func- 
tion of mass is significantly suppressed in all of our 
simulation sets (see Figures ErE]). Summarizing Fig- 
ures HHH we find that for a range of halos between 
5 x 10 4 — 5 x 10 5 Mq the suppression of the number den- 
sity of halos at z = 25 is ~ 20% for v\, c = ler v bc, ~ 40% 
for «bc — l-7<7vbo and ~ 55% for Vhc — 3.4er v b c - As time 
goes on, the suppression of the halo mass function due 
to the stream velocity diminishes with a varied rate - for 
example, by z = 15 (the lowest redshift we consider in 
this paper), for Vhc = loVbc the halo mass function ap- 
proaches the case with no stream velocity (i.e., no sup- 
pression), while in the «bc = 3.4cr v b c case the suppression 
factor remains at ~ 30%. 

Our most significant result is the high fraction of ha- 
los that are almost entirely devoid of gas. We have found 
that for the extreme values of the stream velocity we con- 
sider («bc — 3.4<7 v bc), the gas simply does not accrete on 
the dark matter halos - almost all halos below 10 6 M© 
throughout all of our simulation runs, are empty. For 
more typical values of the stream velocity (wbc = lc v bc) 
most of the empty halos are smaller than about 10 5 M© 
at z < 19. Thus, we expect significant delay in the forma- 
tion o f the first ge n eratio n of galaxies; in our companion 
paper iNaoz et all (pOlTl in prep.) we study the mass 
scales of halos that still retain their gas. 

It was suggested in the literature that gas rich ha- 
los and minihalos may play an important role in cos- 
mic reionization, and that they can p r oduce distinct 
21-cm signatures dKuhlen et al] (12006|); | Shapiro et al. | 
2006) ; INaoz fe Barkanal (|2008t ) but see lFurlanetto fe Ohl 
2006)). Minihalos can block ionizing radiation and in- 



duce an over all delay in the initial p rogress of reion- 
ization (e.g.. IShapiro fc Girouxl Il987t iBarkana fe Loeb 

12004: 



20021: llliev et all 120031: IShapiro et al.l 



llliev et al 



2003 iMcQuinn et al.l 120071) . However, our results here 
(and in our companion paper INaoz et al.ll20TTl m prep. J 
suggest that at high redshift the stream velocity effect 



results in a high abundance of empty halos halos below 
~ 10 5 Mf7> Thus, if r eionization started sufficiently early 
(|Yoshida et al.l [20071 ). in patches of the universe where 
the stream velocity is large, there are fewer gas rich ha- 
los that can absorb ionizing photons. Hence, in these 
patches the delay of the reionization caused by miniha- 
los may be less than in regions that happen to have a 
small value of the stream velocity. The feature is clearly 
quantified by using the gas clumping factor. The vari- 
able measure the dumpiness of baryonic structure out- 
side of bound objects. We found that it is very sensitive 
to stream velocity effect (see Figured]). Therefore, not 
only the formation of the first generations of galaxies 
may be affected by the stream velocity effect, but also 
the reionization process may proceed differently in re- 
gions with very different stream velocities . 

We note that recently iBittner fe Loebl (|2011[ ) showed 
that if the ionizing sources were primarily in halos dom- 
inated by atomic cooling (corresponding to > 10 8 M©), 
then the effect on reionization is negligible. However, for 
sources at high redshift, in halos dominated by molecular 
hydrogen the effect on the evolution of the ionizing front 
may be significant. 
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